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We study correlations on the euclidean spacetimes generated in Monte Carlo simulations of the model. In 
the elongated phase, curvature correlations appear to fall off like a fractional power. Near the transition to the 
crumpled phase this power is consistent with 4. We also present improved data of our computations of the binding 
energy of test particles. 



1. Introduction 

The canonical partition function of the dynam- 
ical triangulation model is defined as a sum over 
simplicial manifolds consisting of equilateral four- 
simplices glued together according to triangula- 
tions T, 

Z=Y^ exp(«2iV2). (1) 

r(W4) 

Here A^2 is the number of triangles in the simpli- 
cial manifold consisting of iV4 four-simplices and 
the topology is chosen to be that of . (For 
more details see [|l|.) We recall that the system 
has two phases, a crumpled phase at low K2 and 
an elongated phase at high ki. 

The model is supposed to represent the quan- 
tum gravitational path integral over euclidean 
spacetimes weighted with the Regge-Einstein ac- 
tion, with volume V oc 7V4 and bare Newton con- 
stant Go oc K^""^. As such it should be able to re- 
produce semiclassical Einstein gravity as an effec- 
tive theory. This can presumably be investigated 
using coordinate invariant correlation functions. 
Another test is to see if scalar test particles form 
bound states, with appropriate binding energies 
under nonrelativistic conditions. 

In 1^ we reported on such preliminary bind- 
ing energy calculations. More insight was needed 
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at that time about the spacetimes described by 
the model. For a semiclassical interpretation in 
terms of an effective action it is obviously impor- 
tant if there is in some sense a background space- 
time which resembles the classical . To address 
this question we studied the number of sim- 
plices at geodesic distance r, N'{r). We found 
encouraging scaling of this function in a region 
around the transition between the crumpled and 
elongated phase. (For a different scaling analy- 
sis see 1^.) At the transition we observed that 
N'{r) behaves approximately like a four-sphere, 
oc sin^[(7r/2)r/rm], for distances r outside a short 
distance 'planckian regime' up to intermediate 
distances of order r„i. For larger distances fluctu- 
ations distort the shape of N'{r). This suggests 
that we focus our attention on this intermediate 
distance regime for a possible semiclassical inter- 
pretation of the data. 

2. Correlations 

Let Tlx denote the number of four-simplices 
winding around a triangle x (often called a 
'hinge'). We measured correlations js] 

Co-(00)(r)-[(01)(r)]^ (2) 

where the brackets denote the average over sim- 
plicial configurations, 

(00)(r) = (5f^^i^^), (3) 
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Figure 1. (n^i n-^){d) and {n,-^ l)^(rf). 



Figure 2. CR{d) corresponding to fig. 



with Ox 



and similar for (01). Here 



dxy is a 'triangle geodesic distance' between x and 
y. The 'disconnected part' subtracted in has 
the form of a correlation of the observable Ox 
with the distance to a point y, dxy, which is itself 
a function of the geometry. This at first sight 
unfamiliar form is needed for the correlation to 
vanish at larger distances. 

Other expressions for similar correlations are 
conceivable. For example, including a volume el- 
ement cx Tlx with a summation over x may ap- 
pear more natural from the Regge calculus point 
of view. We shall argue below that this leads to 
the same scaling results. Taking the average ( ) 
separately in numerator and denominator is an- 
other option, which should not make significant 
difference for not too large distances, because of 
the self-averaging effect of the summations in (||) . 

In Regge calculus Ux^ is related to the scalar 
curvature Rx, and Ux to a local volume element 
at X] we write therefore Cr (Cy) for Ox = 
{Ox — Tlx). However, generally, in a scaling re- 
gion, a lattice operator is expected to reduce to 
a combination of various continuum operators, 
weighted with powers of the lattice distance ac- 
cording to their dimensions. Assuming that the 
theory can be described by a continuum metric 
tensor with corresponding curvature R, the 
dominant continuum observable corresponding to 



both our lattice correlation functions Cr and Cv 
would be given by 



(/ dx dyyM^y/^Sjdix, y) - r)R{x)R{y)) 
(/ dx dyy/g{x)y/^S{d{x,y) - r)) 



{J dxdyy/g{x)y^^S{d{x,y) - r)R{x)) 
(/ dx dy^/g(x)^/g(jj)6{d{x, y) - r)) 



Here d{x, y) is the geodesic distance between the 
points X and y for a given metric g^i/, and the 
average over geometries is supposed to be calcu- 
lable in semiclassical fashion. Of course, we do 
not know the effective action specifying this aver- 
age. It could be a combination of J dx^R and 
higher order terms like R^ . 

Fig. |l| shows the two terms in eq. (||) separately 
for Cfl(r) for Ni = 16000 and K2 = 0.80 (crum- 
pled phase), 1.22 (transition) and 1.50 (elongated 
phase). The two contributions are very similar 
and fig. ^ shows their difference, Cfl(r), which 
drops to zero before finite size effects take over. 

Fig. I shows CrXt) for a system of = 32000 
simplices and K2 = 1.255, which is at the transi- 
tion between crumpled and elongated phase. For 
this case, with the employed definition of geodesic 
distance the 'planckian regime' mentioned in 
the Introduction is r ^ 7, while « 11. So 
the approximate S"* background geometry is to 
be found in the region 7 < r < 11. 
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Figure 3. Power fit to Cfl(d) near the transition. 



Figure 4. CR{d) in the elongated phase. 



We see an approximate power behavior Cy ~ 
ar^ for r > 7. A fit in the range 8 < r < 18 
gives a = -0.5(2), b = -4.0(2). The plot for 
Cv looks similar and leads to a compatible power 
b = -4.30(12), with a = -5.7(1.6). The equality 
of the powers supports the argument given above 
that Cr and Cy are dominated by the same op- 
erators in the scaling region. 

The power law behavior of the data may be 
expected for a flat spacetime background, with 
Cyfl cx Gr~* or cx G^r~^ for one- or two-graviton 
exange. The approximate interpretation of the 
spacetimes near the transition implies finite size 
corrections to the power law. An background 
curvature R would also allow a behavior G^Rr"'^, 
by dimensional reasoning. 

We measured Cv,r also in the elongated phase 
at K2 = 1.50, where the system is known to be- 
have like a branched polymer with effective di- 
mension 2. As fig. ^ shows, a power law looks bet- 
ter in this phase and a fit in the region 3 < r < 15 
leads to 6 = -2.56(3) for C^^ and b = -2.57(2) for 
Cy. Perhaps this can be interpreted in an effec- 
tive conformal mode theory In the crumpled 
phase the correlation data are difficult to anal- 
yse because the distances are rather short in this 
phase. 

We draw the optimistic conclusion from the 
power behavior that the correlation functions in- 



dicate the presence of massless excitations. 

3. Binding 

The behaviour of test particles in the space- 
times generated with the dynamical triangulation 
model is another interesting subject. Test parti- 
cles have no backreaction on the geometry and 
a calculation of their attraction due to fiuctua- 
tions is analogous to the quenched or valence ap- 
proximation in lattice QCD. Let Gxy be a scalar 
field propagator on a given simplicial manifold, 
(-□ 



xy 



Sxz, where now x,y de- 
note the points on the dual lattice and □ is the 
corresponding lattice laplacian. We measured 



G(r) 



G(2)(r 



, ^y ^d^y-r Gxy ^ 

X/j^ ^d^y-r Gxy 
^y ^d^y-r 



(4) 



where dxy is the (dual lattice) geodesic distance 
between four-simplices. On a fiat background one 
would expect the behavior G{r) cx r"exp(— mr), 
G'*^^-'(r) cx r'^exp(— Mr) at large distances, with 
m the (renormalized) mass of the scalar particle 
and M the mass of a two-particle state (in the 
quenched approximation). This is a bound state 
if M < 2m, with binding energy Eb = 2m — M > 
0. 
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Figures. G(2)(d) and G(d)2 



Figure 6. Effective binding energies. 



Fig. I sfiows G(r)2 and G'-^\r) for = 32000 
at the transition, k,2 — 1.255, for several values of 
the bare mass, mo — 1, 0.316, 0.1, 0.0316. The 
approximate S*^ behavior mentioned in the Intro- 
duction is now found in the region 10 < r < 21, 
with the distance definition emplyed here. We 
see roughly exponential behavior suggesting that 
the character of the background geometry can 
in a first approximation be neglected, similar to 
the correlation functions above. The correspond- 
ing renormalized masses are m ~ 1.19, 0.57, 0.27, 
0.14. The ratio m/mo appears to increase sub- 
stantiaUy as mg gets smaller. 

If the model describes semiclassical gravity we 
may expect hydrogen atom like behavior in a 
nonrelativistic situation, with small fine structure 
constant a — > Gw?: 
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4 



(5) 



In this way we may be able to define the renormal- 
ized Newton constant G. To get a rudimentary 
feeling for the corrections to (|^) the hamiltonian 
H = + — Gw? /r may be useful. Re- 

placing p —f 1/r and minimizing the energy leads 
to Ef, = 2m — 2m-\/l — G'^m'^/A, which suggests 
that Gm^ = 2 has to be considered 'large'. 
Fig. H shows the effective binding energy 

off 1 G(r^2 



E[ 



r GC^){r)' 



(6) 



Clearly, there is no indication of nonrelativistic 
rrt' behavior. To judge the situation we have to 
make sure that Gm? is not too large (assuming G 
exists in the model) , or m too small for neglecting 
finite size S"' corrections, or both. It is also de- 
sirable to get a semiclassical understanding of the 
quenched approximation. Finally, we may have to 
take into consideration that the euclidean effec- 
tive action at intermediate distances (neglecting 
higher order terms like E?) is unbounded. 
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